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MÉMOIRE 


Sur le calcul numérique des Intégrales définies . 

Par M. POISSON. 

Lu à l’Académie des Sciences, le 11 décembre 1826. 


Le calcul des intégrales définies est peut-être la partie 
de l’analyse dont les applications sont les plus nombreuses 
et les plus variées. Non-seulement elles comprennent la 
rectification des courbes, l’évaluation des surfaces et des 
solides, la détermination des centres de gravité, mais encore, 
la plupart des problèmes de mécanique ou de physique qui 
se résolvent par le calcul intégral ^conduisent à des expres¬ 
sions des inconnues en intégrales définies. Aussi, depuis 
liuler et surtout dans ces derniers temps, les géomètres se 
sont-ils beaucoup occupés d’étendre et de perfectionner cet 
important calcul. Dans le petit nombre de cas où l’intégrale 
générale est connue sous forme finie, on en déduit immé¬ 
diatement lintégrale définie; dans d autres cas, beaucoup 
plus étendus, on parvient h trouver la valeur exacte de l une 
sans connaître celle de 1 autre; mais le plus souvent on est 
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obligé de recourir aux méthodes d’approximation. Celles-ci 
consistent en des moyens particuliers à quelques intégrales, 
d’après lesquels on parvient à les faire dépendre les unes 
des autres et à les réduire en tables, ainsi que M. Legendre 
l’a pratiqué à l’égard des transcendantes elliptiques et de deux 
autres classes d’intégrales que notre confrère a nommées 
Eulériennes. Quelquefois aussi, on peut réduire la quantité 
soumise à l’intégration, en une série convergente dont les 
termes sont intégrables par les règles ordinaires. Mais quand 
toutes ces ressources manquent, on emploie un procédé gé¬ 
néral de calcul, fondé sur la nature même des intégrales, 
et que l’on appelle proprement la méthode des quadratures ; 
dénomination qui lui vient de ce que le problème est le 
même que celui de trouver l’aire d’une courbe plane, ou le 
côté du carré équivalent. C’est cette méthode, envisagée sous 
un nouveau point de vue, qui est l’objet principal de ce 
Mémoire. 


(i) Une intégrale définie est la somme des valeurs de la 
différentielle, comprises entre les liniites de l’intégration, et 
supposées toutes infinin^it petites, ce qui ne souffre d’ex¬ 
ception que quand le ctjsfficient différentiel devient infini 
entre ces limites. II en résulte que si l’on prend seulement 
un grand nombre de ces valeurs, et qu’on y remplace la 
différentielle de la variable par sa différence finie, on aura 
une valeur de l’intégrale, d’autant plus approchée que cette 
différence sera plus petite; et la méthode dont nous allons 
nous occuper, consiste à trouver, exactement ou par ap¬ 
proximation, la correction qu’il faudra faire subir à ce pre¬ 
mier résultat. 
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Nous supposerons les valeurs de la variable qui repondent 
aux deux limites de l’intégrale, égales et de signes contraires, 
ce qu’on peut toujours obtenir en augmentant ou diminuant 
la variable d’une quantité constante. Nous désignerons ces 
limites par =b a; et pour les indiquer en même temps que 
l’intégrale, nous emploierons la notation très-commode que 
M. Fourrier a proposée. Ainsi 



désignera l’intégrale de fxdx, prise depuis x =— a jus¬ 
qu’à x= a ; f x étant une fonction donnée qui ne devient 
pas infinie entre ces limites. 

Partageons a en un nombre n de parties égales; soit w la 
grandeur de chacune d’elles,en sorte qu’on ait a= Ru ; fai¬ 
sons, pour abréger, 

"./*( 71 <*> ) + — 71 co -f- Ci) ) Ci) -f- 2 Ci) ) -+-. 

. f ( TL CD -- 2(i>) -J- Tl a)- ti) ) -t-/(n M ) = P„: 

en remplaçant dx par o, on pourra prendre, d’après le prin¬ 
cipe précédent, wP„ pour la valeur approchée de notre inté¬ 
grale; et si l’on désigne par Q„ la correction dont elle est 
susceptible, on aura exactement 


r a 

I fxdX — wP« + Q«. 


(0 


Au lieu de ne faire entrer dans P„, qu’une seule des deux 
valeurs extrêmes de fx, on a pris, pour la symétrie du cal¬ 
cul , la moitié de chacune d’elles ; ce qui est permis tant que 
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la correction Q* n’est pas déterminée. Il s’agit maintenant 
de trouver la valeur inconnue de cette quantité, en fonction 
du nombre n , ou de la différence w. 


(2) Pour y parvenir, je ferai usage de la formule : 

f x =i~ a ,f f x ' dx '+ \ /_Js7 cos. dx', ( 2 ) 


qui se trouve dans mon dernier Mémoire sur les intégrales 
définies, inséré dans le dix-neuvième cahier du journal de 
l’Ecole Polytechnique. Elle représente fx pour toutes les va¬ 
leurs de Æ‘> — a et <a; à chacune des valeurs extrêmes 
x—zka, le second membre de cette équation est égal à la 
demi-somme des valeurs correspondantes dey#, c’est-à- 
dire, qu’il faut joindre à lequation (2), celle-ci : 


[/«*+/(—' a )]=-taf f x ' dx ' +l a[ [2 cos - '~ (a fl X) \fx dx'. 
J —a J —a>- 1 


On représente, à l’ordinaire, par r le rapport de la circonférence 
au diamètre; i est un nombre entier et positif, et la somme 

comprise sous les signes J', s’étend, comme l’indique la no- 

00 

tation 2 , à toutes les valeurs de i depuis i— 1 jusqu a 

1 

1 = oc . Mais avant d’employer ces équations, il convient de 
rappeler la démonstration que j’en ai donnée, en regardant 
chacune des séries périodiques quelles renferment, comme 
la limite d’une autre série convergente. 

Soit donc a une quantité positive et moindre que l’unité. 
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Considérons l’expression 


7 


S x ' dx ' + 'af [2 “'COS. ~\fx dx'; 

J —a J -a'- i J 


le second membre de 1 équation (2) sera la valeur de X qui 
répond à la limite où la différence 1 — a est infiniment pe¬ 
tite ; ainsi il s’agit de faire voir qu’à cette limite on a X —fx, 
pour x > —a et < a . 

Or, en développant suivant les puissances de a, on a, en 
série convergente, 


- 7 . 0 . cos. —-- + a’ 


= 1+22 a'COS. 


ir:(x -~æ J ) . 


on aura donc 


r_ i Ç a ( i — o'^fx'dx' 

L 2 a I — x') 

J —a 1 — 2 a cos. —- ’ 


+ a* 


Le coefficient de dx' sous le signe intégral, devient infini¬ 
ment petit en même temps que 1 — a , excepté pour les va¬ 
leurs de x qui rendent cos. — ■ ~~* ^ infiniment peu diffé¬ 
rent de l’unité, et, par conséquent, le dénominateur aussi 
infiniment petit. Mais x étant > —<2 et < a, et la variable 
x ' ne sortant pas non plus de ces limites, cette circonstance 
ne peutavoir lieu que pour des valeurs de x '— x, infiniment 
petites, positives ou négatives; il suffira donc d’étendre l’in¬ 
tégrale aux valeurs de x' infiniment peu différentes de x; et 
dans cette étendue, on pourra considérer la fonction fx 
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comme constante et égale à fx . Si donc on fait 


i — *=g, 


x —— x "H A ^ d x — d h f 


et que l’on traite g et h comme des quantités infiniment pe¬ 
tites , on aura pour la limite demandée : 


X 


/f Ç § dh 
' aJ g •+=^‘ 


Comme cette dernière intégrale est infiniment petite, pour 
toute valeur finie de la variable, on pourra l’étendre à des 
valeurs de h aussi grandes que l’on voudra ; en désignant 
donc par £ une quantité positive et finie, dont la grandeur 
est arbitraire, et intégrant depuis h = —£ jusqu à h= -+- &, 
nous aurons 


X=^arc.(tang. = |); 


quantité qui se réduit à X=fx, à cause de g infiniment 
petit ; ce qu’il s’agissait de démontrer. 

Dans le cas de x=a, on rendra cos. 77 ^ g X ^ infiniment 

peu différent de l’unité , en supposant successivement 
x'=a + h et x' = —a -\- A, et traitant toujours h comme 
une variable infiniment petite; mais pour que x ne sorte pas 
des limites ±a, il faudra n’intégrer que depuis h — — $ 
jusqu’à h = o dans la première hypothèse, et depuis h= o 
jusqu’à h= + £ dans la seconde ; ce qui réduira chaque por¬ 
tion d’intégrale à la moitié de la valeur précédente. Alors 
on aura, dans ce cas, 
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pour la limite de X ; et l’on trouvera le même résultat dans 
le cas de.z = — a; ce qui coïncide avec l’équation ( 3 ). 

( 3 ) Maintenant, pour faire de l’équation (a) l’usage que 
nous avons en vue, mettons-y nuà la placé de a \ puis don¬ 
nons successivement h x, ces 2 n —1 valeurs équi-diffé- 
rentes : 

— {n — i)w,^-(rc — 2)w... —to,o, -h «,. .( n —2)(û,(n— i)w, 

poun lesquelles cette équation subsiste. En prenant la somme 
des résultats, et ajoutant la demi-somme des valeurs defx 
relatives k n = =t a, on en conclura 

» p .= ;[/ Æ +/(- a )] + 2 " ~ ‘ f fx'dx 

' —a 

+ "/-« P COS ’ 4 ? ]/*’**' » 

où l’on a fait, pour abréger, 

ÏTC ’li'K 3i7T (n _iV-tt 

p= I + 2 COS. h 2 COS.—< h 2 COS.-h . . -h 2COS.^- -2—- 

1 n n n n 

Cette quantité p est évidemment égale à in — i toutes les 
fois que? est un multiple de 2 n. Si on la multiplie par cos. — , 
on trouve, en réduisant, 

* tc . (n — l'l i 'K 

p cos. — =p 4- cos. I ir — cos. i- - — ; 

d’où l’on tire p= —cos.jtt, pour les autres valeurs de i. 
D a près cela, je prends d abord la somme 2 en supposant 


2 
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que cette dernière valeur de p ait lieu sans exception; en¬ 
suite, je fais i = 2Ïn, et je prends une seconde somme 

00 

2 relative à i\ en supposant que 2n — i 4- cos.2«'«7r, ou 

2/î, soit la valeur de p : la série périodique, contenue sous 
le signe intégral, se composera évidemment de ces deux 
séries ainsi calculées ; et en faisant attention qu’on a 


cos. I tz cos. 


*7ç(a — x') 


on en conclura 


v ir x’ 

2 P cos ——=—2 cos. 

Tl ft) 


iizfa — x') 00 li'izx' 

—-- -4-2 n\ cos.- 


De cette manière, la valeur précédente de wP„ deviendra 


»P.=ll/a+/(—a)l+ (i—l)j a fx’dx' 

En ayant égard à l’équation ( 3 ), celle-ci se réduit à 
»\>.=f a _ fxdx +zf a [ 2 *C 08 .Î^] fx'dx'-. 


et si 1 on y remplace x et i par x et i, et qu’on la com¬ 
pare à l'équation (i),il en résultera 
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Ainsi la correction quon doit faire subir à la première 
valeur approchée de notre intégrale, se trouve exprimée par 
une autre intégrale définie; mais, par le procédé de l’inté¬ 
gration par partie, celle-ci se réduit en une série ordonnée 
suivant les puissances de <o, dont il suffira généralement de 
considérer les premiers termes. 

( 4 ) Soit, pour cela , 


i +^ + 3T^ + 4T7 n -+- et c*=T(a7r) ll "A m ; 

la série se prolongeant à l’infini, m étant un nombre entier, 
et A „ un coefficient dont la valeur exacte est connue pour 
toutes les valeurs de m. En intégrant 2m—i fois de suite 
par partie, et observant qu’aux limites dza, ou ±ft W , on a 


HT X • HT. X 

cos. -=i, sm. --= o. 


nous aurons 


Q:=- 




\ L dx J V. dx t 

-•([î£H3£)>- 




les différentielles comprises entre les parenthèses se rappor¬ 
tent à la première limite x = — a, celles qui sont renfer¬ 
mées entre des crochets répondent à la seconde x=za, et 


2. 
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l’on a fait, pour abréger. 


».—<-,)• (ÛÏ'C ££<<* ( 4 ) 


Maintenant, si l’on prolonge indéfiniment cette série , et 
qu’on néglige le reste R„, l’équation (i) deviendra 

f fxdx 

J — a 


Cette formule ne différé pas essentiellement de celle 
qu’Euler a donnée pour le même objet dans son Traité de 
calcul différentiel, et qu’on peut regarder commeune des plus 
utiles dont il a enrichi l’analyse; mais la manière dont nous 
y sommes parvenus a l’avantage de faire connaître en même 
temps l’expression du reste R m , que l’on néglige quand on 
s’arrête à un terme quelconque de la série infinie ; expres¬ 
sion dont il sera toujours facile d’assigner une limite supé¬ 
rieure à sa valeur exacte ; ce qui permettra d’apprécier le 
degré de l’approximation. Il serait à désirer que l’on efrt de 
semblables limites pour toutes tes suites infinies dont on fait 
usage : Lagrange les a exprimées très-simplement dans le 
cas du théorème de Taylor; et récemment M. Laplace s’est 
occupé de questions analogues, relatives aux développe¬ 
ments des coordonnées des planètes dans le mouvement 
elliptique, et d’une autre fonction qui se présente dans la 
théorie des perturbations. 


—p—■([&]-(£)>. ' 


\ L dx 
etc. 


( 5 ) 
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( 5 ) Les valeurs des coëffîcients A,, A a ,A 3 , etc., sont con¬ 
nues, comme nous l’avons dit ; mais on peut aussi les déter- 
miner au moyen de l’équation ( 5 ), en faisant une supposi¬ 
tion sur le nombre n et sur la fonction fx. Le plus simple 
est de prendre n— i. On a alors 

<a=a, P„ = 7 t /(— a) 4 - fo 4 - '-fa. 

Si l’on désigne par e la base des logarithmes népériens, et 
qu’on fasse 

fx = e*, 

on aura 




et généralement 


r d im -'fx i_ (_ 

L dx 2m -' J \dx' m - 1 ) 


a 


e —e 


au moyen de quoi, en divisant les deux membres de l’équa¬ 
tion ( 5 ) par e a — e ~ a , on en conclura 


a (^c ' a —e 1 


-h a' A I — a* A, -t- a 6 A 3 — etc. 


Or, le premier membre ne changeant pas quand on y change 
le signe de a, son développement ne renfermera que des 
puissances paires de a; d’ailleurs, les inconnues A,, A„ A 3 , 
etc., sont indépendantes de a ; il en résulte donc qu’après 
avoir développé le premier membre, les coëfficients de 
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a 3 ,— a\ etc., seront les valeurs de ces quantités; ce 
qui était déjà connu. 

En supposant toujours n = i, (o=a f et prenant 
on aura, m étant un nombre entier, 


/ a 

fxdx = 

-a 


2 m -\-1 


d’après l’équation ( 4 ), le reste R m sera nui ; et en supprimant 
le facteur qui se trouvera commun à tous les termes 

de l’équation ( 5 ), il vient 

— ■ — a m A, — 0 . 771.1771 — 1.2 772 — 2 . A a -+-. 

2 (a m -+- i) 

.± 2 m . 2 m — 1.2 m — 2. . . 3 .2. A,„. 


Cette relation connue entre les quantités A,, A,, A 3 , etc., 
servira à les déterminer les unes au moyen des autres, en y 
faisant successivement m= 1 , = 2 , = 3 , etc. On trouve, de 
cette manière* 





etc. 


( 6 ) Si l’on veut que l’intégrale proposée commence avec la 
variable, et soit prise depuis zéro jusqu a une limite donnée c, 
on diminuera x de a, et l’on fera 2Mettons ensuite 
y#à la place de f(x —<2), et n au lieu de 2 n; la différence 
w sera la n. ième partie de c ; la quantité P, se composera de 
n -h 1 termes, savoir : 

P„= i/o +/c> 4-/2w + . . . +/(«<!> — a) + ifc; 
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et si l’on substitue à la place de A t , A,, A 3 ’, etc., leurs Valeurs 
numériques, l’équation (5) deviendra 

_ ^ (\ d ü^~\ ( d 'f x X\ 

3oa4o\L^j: 5 J \dx* JJ 

-h etc. ; 

les différentielles comprises entre des parenthèses ou entre 
des crochets, répondant toujours, les premières à la pre¬ 
mière limite x — o , et les secondes à la seconde limite x=c. 

A cause que a est un multiple de w, le cosinus de 
reste le même dans le changement de x eu x — a; le reste 
donné per l’équation (4), qu’il faut ajouter à cette série 
quand on s’arrête au n. ième terme inclusivement, aura donc 
pour expression : 



R *=~»(-*)'(£)*"/ h cos -^ll^ dx :S?) 


En intégrant encore une fois, et observant que sin. 


est nul à la seconde limite x = c=n<ù , aussi bien qu’à la 
première x = o, on peut donner à R,« cette autre forme équi¬ 
valente : 





dx 

dx'»+' 


• (8) 
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Pour obtenir des limites de ces expressions, désignons 
par B„, et C m des quantités connues, telles que l’on ait, abs¬ 
traction faite du signe y 


d* m fx 
dx lm 




fx 
d 


< 


C m , 


pour toutes les valeurs de x comprises depuis zéro jusqu’à 
c; soit, en outre, 


, / \ 2W 4" I a / 

etc. =7(27:) • A 


en observant qu’on a évidemmment 



d’après la formule (8). * 

Ces limites montrent que pour une même valeur de m , 
ou pour un même nombre de termes de la série (6j, l’ap¬ 
proximation croîtra indéfiniment à mesure que w diminuera ; 
mais quel que soit <*>, l’approximation ne croîtra pas de même 


P 
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avec le nombre m, et il arrivera très-souvent qu elle décroîtra 
au-delà d’un certain nombre de termes. C’est ce qui aura 
lieu quand les quantités B,„ et C w augmenteront plus rapi¬ 
dement avec m, que &>*'" ne diminuera, et alors la série (6), 
après avoir été convergente dans les premiers termes, de¬ 
viendra divergente et par conséquent inexacte. Dans le cas 
de c = oo ,ces limites seront illusoires, et il en faudra déter¬ 
miner d’autres , propres à chaque exemple particulier. 

La valeur exacte de la quantité A' m qui entre dans la se¬ 
conde limite, n’est pas connue comme celle de A,*. On a 
fait usage de différentes méthodes pour calculer sa valeur 
approchée; nous indiquerons celle que fournit l’équation(6), 
en y faisant 

fx — -,—-— \ t , «=i, c=oo. 

J (i-h x) ‘ ’ - 9 

On a alors 


J fxdx= p» = — 7+ ;(2*) 1 " i + 1 A' m ; 
et quel que soit le nombre entier i, on a aussi 

Au moyen de ces valeurs, on tire de l’équation (6) : 


( aK ) 2 ' ra+, A' 


m m 


w + i t m r m H- r. m -4- 2 . m -h 3 

~6 36*o 


n. -|- 1 . m + 2 . m -4- 3. m -f- 4. m 4- 5 
1 5 1 2 o 


• etc. 


En appelant y. le nombre de termes de cette série que l’on 
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conservera, à partir du second, et M le reste qu’il y faudrait 
ajouter, on aura, d'après la formule (7) : 


1 ) fA + 1 ,m + 2.. .m+ ap.J' [2 j7jzC08.2ïwa?|j| 

Il sera alternativement positif et négatif, ce qui montre que 
la série précédente donnera des valeurs de A f m , alternative¬ 
ment plus grandes et plus petites que la valeur exacte, et 
qui en différeront, par conséquent, d’une quantité moindre 
que le terme où l’on s’arrêtera. 

(7) On peut éliminer les différentielles de la fonctiony# 
qui sont contenues dans la formule (6), et les remplacer par 
ses différences finies. 

En effet, soit 

Fz=/z +f(c — z) ; 

désignons par A 0 A,, A 3 , etc., les différences successives de 
F z, qui répondent à z = o, et sont prises en supposant 
Az = u, de sorte qu’on ait 

a,=F 2« — Fo, 

A^Fûw — 2Fw + Fo, 
a 3 =F 3 w— 3 F 2 w -h 3 F (ù —Fo, 
etc. ; 

nous aurons cette formule d’interpolation : 


Fz = F o -4- - A, 


z(z — (ù) (z- 2 <*>) 


d’après les notations précédentes, on aura aussi 
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fd' m ~'fx\ f d""-'fx\ d' n - x Y z 

Kdcc"»-' ) L^ — “J- d*'—' ’ 

pourvu qu’on fasse z= o après la différentiation : il en ré¬ 
sultera 




etc. ; 


et en substituant ces valeurs dans l'équation (6), il vient 


£ fx dx — wP n -+- A 


G) . IQb) 

—7 A, H- 

24 720 


A 3 


3 Gi 

t6o 


A, + etc. (9). 


En y faisant g>= 1 , cette formule coïncide avec celle que 
M. Laplace a donnée dans le tome IV de la Mécanique céleste. 
Elle ne suppose pas connue l’expression de fx. Pour en faire 
usage, il suffira d’avoir un nombre n H- 1 de valeurs numé¬ 
riques de cette quantité, correspondantes à autant de valeurs 
équi-différentes de x: on ne pourra toutefois l’employer uti¬ 
lement que quand les différences A,, A a , A 3 , etc., décroîtront 
très-rapidement. 

La formule d’interpolation dont nous sommes partis, ne 
subsiste pas lorsque/# est une fonction périodique desin. 27c x 
et cos.27r x, et que l’on prend w= 1 ; car alors toutes les 
différences A,, A,, A 3 , etc., quelle contient, seraient milles, 
et l’on aurait Fjz=Fo, ce qui n’est pas vrai. Il en résulte 
que l’équation (9) que nous en avons déduite, n’aura pas lieu 
non plus dans ce cas particulier; mais on obvie à cet incon- 

3. 
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vénient en prenant une autre valeur pour la différence w qui 
peut être aussi petite que l’on voudra. 

(8) La formule d’Euler se trouve en défaut dans un autre 
cas, ainsi que M. Legendre en a fait la remarque à la fin de 
son traité des Fonctions elliptiques (i). Ce cas a lieu lorsque 
les différentielles impaires d efx s évanouissent, ou, plus 
généralement, sont égalés aux deux limites de 1 intégrale 
que l’on considère. Quelle que soit la quantité <*>, 1 équa¬ 
tion (6) se réduirait alors à 


y fxdx-=toV n ; 

O 


d’où il résulterait que l’intégrale proposée s’exprimerait sous 
forme finie, et que sa valeur dépendrait du nombre arbi¬ 
traire n, ce qui serait absurde. Mais ce résultat provient de 
ce que l’on a négligé le reste Rm, qui, dans ce cas particu¬ 
lier, au lieu de décroître h mesure que m augmente, est au 
contraire indépendant de la grandeur de ce nombre. En 
effet, en intégrant par partie, et observant qu’on a, par 
hypothèse, 


_ [d' m 

d.r 2m ~ hI J L dx 



la formule (8) donnera 

jn/üNam + i/ rr” i _ a i Ttx'] d' m + \f x A 

Uï) J Q 11, î^ cos -—J m^ dx 

et si l’on compare cette expression à la formule (7) dans la- 


*(r) ; Tome II, page 578. 
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quelle on augmentera m d’une unité, on en conclura 

R =R ; v 

c’est-à-dire, que le reste R™ est constant par rapport au 
nombre m. 

Dans ce cas singulier, si ^ est une fraction, le facteur 

m °I ue ce reste ren f erme i diminue indéfiniment à 

mesure que m augmente; mais en même temps l’autre fac¬ 
teur augmente à raison des différentiations successives de 
fx , de telle sorte que le produit demeure constant. Cepen¬ 
dant, il n’en sera pas de même à l’égard de la limite supé¬ 
rieure à R,„ que l’on pourra assigner; elle dépendra de m, 
et décroîtra d’autant plus rapidement avec w que ce nombre 
m sera plus grand. 

(9) Nous choisirons pour exemple de cette anomalie, l’in¬ 
tégrale qui donne le quart de la circonférence d’une ellipse, 
que nous appellerons E. En prenant pour unité le demi- 
grand axe, et désignant l'excentricité par h x nous aurons 
alors 

E = / t/1 — h * siu. a x dx ; 


on fera donc, dans l’équation (6), 

£=777, f x-=± \/x— iSsin.'x; 

et il est évident que toutes les différentielles impaires de fx 
seront nulles aux deux limites x = o et #=7*, à cause du 
facteur sin. x cos. x dont elles seront affectées. 
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Le coefficient différentiel de fx, d’un ordre quelconque, 
se composera de termes dont les diviseurs seront les puissan¬ 
ces impaires de fx, et qui ne renfermeront quedes puissances 
entières de sin. 2# et cos. 2 x h leurs numérateurs ;-on substi¬ 
tuera facilement à chaque numérateur, une plus grande quan¬ 
tité abstraction faite du signe ; prenant en outre, pour fx sa 
moindre valeur \/1 —À*, on formera ainsi une quantité su- 

périeure au coefficient différentiel ^ x 7 ^r entre dans la 

formule (7). Soit H cette quantité; en mettant aussi dans 
cette formule , 7(2*)»^,* à la place de la série quelle ren¬ 
ferme , et pour (o sa valeur ^, nous en conclurons 



en grandeur absolue. 

Appelons I le rayon du cercle dont la circonférence est 
équivalente à celle de l’ellipse que nous considérons, de sorte 
qu’on ait E=7xl. D’après la formule d’Euler, nous aurons 
cette valeur approchée : 


I=- ( i/o +/—H-/— -+-/—■+..+/- 

«V * 7 J in J 2 « J 2 n '' 


À n — 0 * 


;/*) ! 


et si l’on désigne par $1, l’erreur dont elle susceptible, on 
aura en même temps 


«<(£)*■ 


HA 


En supposant ^=o,6, et s’arrêtant à m = 3 , on trouve 
qu’on peut prendre pour H, le nombre 4 °> et ^ on s ub- 
stitue, en outre, pour x et A 3 leurs valeurs numériques, il 


1 
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vient 

H ■■<&£&. 

Soit ensuite n = 4 ; la valeur approchée de I sera 

i = o,9Q 2 7799 2 7 2 ‘> # 

à quoi il faudra joindre 

£1< o,ooooo 483 i i; 

de sorte que l’erreur sera moindre qu une demi-unité déci¬ 
male du 6 e ordre. Elle est, en effet, plus petite; car en 
calculant, par les méthodes particulières aux transcendantes 
elliptiques, une valeur de I exacte jusqu’aux décimales du 
i 4 e ordre inclusivement, M. Legendre a trouvé quelle ne 
différait de la précédente que d’une demi-unité du io^ ordre. 

11 est à remarquer que dans ce même exemple, la série 
des différences, contenue dans la formule (9), n’est pas con¬ 
vergente; et qu’en y ayant égard, on s’écarte plus de la va¬ 
leur exacte de E, qu’en s’en tenant à la seule partie wP„ de 
cette formule. 

(10) Si l’on supprime en entier dans l’équation (6), la série 
des différentielles de fx, il faudra ajouter à son second 
membre la valeur de R,* qui répond à m= o. Je fais passer 
cette quantité dans le premier membre, et je suppose qu’on 
ait c = ao, ce qui change uP n en une série infinie; il en 
résulte cette transformation d’une série dans une autre : 

pfxdx+î-Z 1 (f™cos.î^fxdx')=-,fo+ ^/‘“I ( 10 ) 
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en déduisant l’expression de R G , de la formule (7). Mais pour 
que cette nouvelle formule soit utile, il faudra l’appliquer 
à des cas dans lesquels les intégrations relatives à x puissent 
s’effectuer sous forme finie. 

Si nous prenons, par exemple, 

*fx-=.e * cos. 2. a x, 

e désignant la base des logarithmes népériens, et a une 
constante donnée, nous aurons 


* cos - v ^rf xdx = ^ [« (" + -) +e ( a .*) j; 

d’où nous conclurons 


l/irze V co/=<o2e cos. 2.ia*>; 

les sommes 2 s’étendant actuellement à toutes les valeurs 
de i, entières, positives, négatives ou zéro, depuis i ——00 
jusqua i= 00. Cette équation est identique dans le cas de 

<z = o et 77. En faisant (o==\y / ^jz 1 a ===-^—= 7 et, pour 

abréger, e ~* = £, on en déduit 

.(H'-n-oV, 

ou, ce qui est la même chose. 


1 i i ü ü 1 

-e 4 4 - 6 * +e < + t‘ -f- etc.=- 1 -h 8 4 —e’-+- etc.: 

2 a 


équation entre les deux transcendantes e et 77, qui mérite 
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d’être remarquée, quoiqu’elle ne soit pas sous forme finie, 
et qu’il ne semble pas qu’on puisse ly ramener. 

(n) Soit encore 

^ cos. a .v 

J X — b '+x' ; 

a et b étant des constantes que nous regarderons comme 
positives. On aura, par les formules connues, 


/. 


00 2 cos. 2 inx cos. a x 




— h(*iiz—a) — a) 


y (») 


le signe supérieur ou le signe inférieur ayant lieu dans la 
première exponentielle, selon que Ion a a < ou > 2 ïtc, afin 
que son exposant soit toujours négatif. D’après cela , si l’on 
désigne par imz, le plus grand multiple de qui soit 
contenu dans a , il faudra prendre le premier signe, lors¬ 
qu'on aura i > n, et le second , dans le cas de i—n ou < n. 


En partageant la somme 2 en deux autres, dont l’une soit 

i 

prise depuis i= i jusqu’à i=n, et l’autre depuis i 

jusqu’à i= oo , et sommant les progressions géométriques 
qui en résulteront, on en conclura 



2 COS. 2 Itz X COS. d X 
h'+x' 


dx= 


— h(2i r+2/zi:— a) — b(ïTZ + a) 

tc e K ; + e 


2 b 


- 2 Tzb 


On a d’ailleurs 


b(zn — a) Æ(air-h 2 « 7 T — a) 
tc e K '—e 


2 b 


2 7C b 



cos .ax i* — 

~b % -\-x* 2 b e 


^ dx 




4 
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si donc on prend = i, la formule (10) deviendra, toute 
réduction faite, 

7 : [ e ~ l '( a - 2n7: -' !Z ) + 2 « 7 T-^ 7 T)j 

777 ^ — rW. 3 \ ~ 

—e J 

I , COS. « COS. 2 Cps.3ût 

ÏP î T +~l + P+4 + 1 F+g + etC ' 

En différentiant cette équation par rapport à a, on en 
déduit cette autre : 



-b{a 


n ( ifb — k£>\ 
2 \e —e ) 


sin.Æ asin. 2 rt 3sin.3a 
/ q — TT~. — T - TT-, - 1 - etc. 



Au moyen du nombre n quelle renferme, lequation (i 3 ) 
subsiste pour toutes les valeurs réelles et positives de a, 
parce qu’en effet, les équations (n) et (12) dont nous som¬ 
mes partis, ont lieu sans exception. Mais il n’en est pas de 
même à 1 egard de leurs différentielles par rapport à a, et 
pour cette raison, l’équation (i4) est en défaut quand a est 
zéro ou un multiple de 9. TC . 

Pour la rendre applicable à ces valeurs particulières, j’ob¬ 
serve qu’en différentiant l’équation (12), et faisant ensuite 
a = o, on aurait 



x sin. ax 
b'+x* 


dx=- 

2 


ir, 


tandis que cette intégrale est évidemment nulle; dans le cas 
de a=o, il faut donc retrancher du premier membre de 
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1 équation (i4), le terme 77: qu’il renferme de trop; ce qui 
le réduit effectivement à zéro, comme te second membre. 
Dans le cas de a = ^n r , et iï= 7 i, la différentielle de l’équa¬ 
tion (11) par rapport à a , donnerait 



2 X COS. 2 

b 


n 7 rxsin. 2 mzx 7 - / 

--3- dx = - i + e 

2 \ 


- 4^ 77 



au lieu que la valeur exacte de cette intégrale est seulement 
\ize—^ n ' r - b \ d’où. il résulte que pour la valeur particulière 
<7. = 2/* 77, le premier membre de l’équation (i4) renferme 
aussi un terme 7 % qui ne devrait pas s’y trouver : en l’en 
retranchant, ce premier membre devient nul en meme temps 
que le second. 

Les équations (ci) et (la), ainsi que leurs différentielles, 
et, par conséquent, les formules (i 3 ) et (* 4 ) qui en dérivent, 
subsistent encore quand on y remplace b par g b\Z~ x ; 
g et b étant des quantités réelles , dont la première est po¬ 
sitive, mais aussi petite que l’on voudra. Après cette sub¬ 
stitution , si l’on suppose que la partie g devienne infiniment 
petite, et qu’on la supprime en conséquence, on aura 


77 cos. b [a — 2 « 77 — 7 ;)_ 1 t cos. a 

ibsxn.Ttb 2 b* b ’—1 

77sin.£(a —2 mz — 77 ) sin.a 

2sin.7r£ b À —i b *—4 


cos. 2 a cos. 3 a 
2sin.a 3sin.3a 


etc., 


etc. 


Toutes les formules de ce n° étaient déjà connues. Elles 
se trouvent dans les ouvrages d’Euler et de M. Legendre, 
et aussi dans mes Mémoires sur les intégrales définies qui 
font partie du Journal de l’Ecole Polytechnique. On en dé¬ 
duit facilement tous les résultats que l’on a trouvés jusqu’à 

4 - 
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présent, et, vraisemblablement, tout ce cjuil est possible 
d’obtenir, sur les séries de sinus et de cosinus, et sur celles 
des puissances négatives des nombres naturels. 

(i2) Soit que l’on forme la valeur exacte d’une intégrale 
définie ou qu’on la calcule par approximation, il faut avoir 
égard aux observations suivantes par lesquelles nous ter¬ 
minerons ce Mémoire. 

i°. Lorsque l’une des limites de l’intégrale est infinie, la 
fonction fx comprise sous le signe j, doit décroître à me¬ 
sure quelle s’en approche, et devenir nulle à cette limite. 
Cela est nécessaire pour que la partie uP« de la formule (6), 
qui se change alors en une suite infinie, soit une série con¬ 
vergente. Néanmoins on a souvent employé des intégrales 
de fonctions périodiques, prises depuis zéro jusqu’à l’infini }> 
mais les valeurs qu’on leur assigne ne sauraient se vérifier 
numériquement, ni être données par la formule (6); et l’on 
doit ne les considérer que comme des limites d’autres inté¬ 
grales pour lesquelles la fonction fx était décroissante et 
convergente vers zéro. C’est ainsi qu’en désignant par a une 
constante réelle et qui ne soit pas nulle, on a 


y cos.axdx = o y f sin .axdx =- 7 
O o 

en regardant ces résultats comme les limites de ceux-ci : 


e êX cos.axdx=- T £—ï, f e ^sin axdx 
1 0 g + a J Q s ~ ha 


dans lesquels est une constante positive, aussi petite qu’on 
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voudra, et que l’on fait infiniment petite pour passer aux 
intégrales précédentes. La première peut encore être consi¬ 
dérée comme la limite de l’une ou 1 autre de celles-ci . 


r 

J O 


_ pr a x* r _ t/ïc 

8 cos .axdx=—e 


f. 


00 


S, 



qui s’accordent aussi à donner zéro pour valeur de cette 
intégrale, à la limite ou 1 on suppose la quantité g infini¬ 
ment petite, en exceptant toujours le cas où l’on aurait 

<3 = 0 . 

A cause de cette exception, si les intégrales que nous citons 
pour exemples, sont comprises sous d’autres signes j, rela¬ 
tifs à a, il faudra avoir soin de les remplacer par celles dont 
elles sont les limites. Ainsi, en désignant par t a> et F a 
deux fonctions données de a, si Ion a 


fFa^f^ïos.axdx^da, JF'a sm.axdx) da, 


1 faudra prendre 


/ gY ad a Çt 

’ J 


'a Y' ad a 




et ne faire la constante g infiniment petite qu’après avoir 
effectué les intégrations relatives à a, lorsque leurs limites 
comprendront a=o. Soit, par exemple , 

Ftf=cos.(3, F'<z=sin.a; 


3o 
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et intégrons depuis a =—co jusqu’à a=oo. Nous aurons 
■ne—s pour la valeur commune aux deux intégrales, laquelle 
se réduira a ir, a la limite ou la quantité g est infiniment 
petite. On peut vérifier que * est, en effet, la véritable va¬ 
leur de chaque intégrale double, en effectuant les intégra¬ 
tions dans un ordre inverse, c’est-à-dire, en commençant 
par a et et finissant par ce qui n’est sujet à aucune dif¬ 
ficulté; ou bien encore, en intégrant d’abord par rapport 
à ar, entre des valeurs indéterminées de cette variable, qu’on 
ne fera infinies qu’après l’intégration relative à a. 

i°. On ne doit pas faire usage de la formule (6), quand 
la fonction^» passe une ou plusieurs fois par l’infini, entre 
les limites de lintégration. Le principe qui en est la base, 
et l’équation (2) dont nous l’avons déduite, supposent essen¬ 
tiellement que y# est toujours une quantité finie. Si cepen¬ 
dant cette fonction devenait infinie à raison d’un diviseur 
dont l’exposant serait moindre que l’unité, il serait facile 
de le faire disparaître par un changement de variable. Sup¬ 
posons, par exemple, 


a étant une constante comprise entre les limites de l’intégra¬ 
tion , k un exposant > o et < 1, et Far une fonction qui ne 
devient pas infinie : on fera alors 

(x — a) l ~ / ‘= r , (1 —k)(x—àf k dx = dy; 

d où l’on conclura 


3i 
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F'y étant aussi une fonction qui ne deviendra pas infinie, 
ce qui permettra d’appliquer la formule (6) à l’intégrale re¬ 
lative à la nouvelle variable jr. 

Généralement, une intégrale f fxdx cesse de représen¬ 
te 

ter la somme des valeurs réelles de la différentielle, com¬ 
prises depuis x=o jusqu’à x — c, lorsque fx devienfrin- 
finie dans cet intervalle. Ainsi, l’on a , par exemple , 

r c dx _| a — c Ç c dx _i i 

J 0 ^ — ° g ' « 7 J Q yx — af~a a — c ’ 

et si l’on suppose a > o et < c, la première valeur est imagi¬ 
naire, et la seconde négative, tandis que la première diffé¬ 
rentielle est toujours réelle, et la seconde toujours positive. 
Mais dans ces sortes de cas, si l’on fait passer la variable, 
entre les limites données, par des valeurs imaginaires qui 
ne rendent plus fx infinie, on pourra considérer de nou¬ 
veau l’intégrale ^ fxdx, commua somme des valeurs ima¬ 
ginaires de fxdx. Dans les exemples précédants, il faudra 
faire 


#=7c(i —cos.z-t-sin.zl/^ 7 ), Ja? = 7c(sin.z h- cos.si/-1), 

et intégrer depuis z — o jusqu’à z=(a/î-h i)^> n étant un 
nombre entier quelconque, afin de ne pas changer les limites 
données x = o et x=^c. Les intégrales definies ne change¬ 
ront pas non plus; mais les fonctions de z, comprises sous 

les signes f, ne devenant pas infinies, chaque intégrale sera 
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maintenant la somme des valeurs de la différentielle; et ces 
valeurs étant imaginaires, cela explique comment leur somme 
peut être imaginaire dans un cas et négative dans l’autre. 

Il y a aussi des intégrales dans lesquelles la fonction Jx 
passe une infinité de fois par l’infini, et qu’on peut encore 
admettre dans l’analyse en les considérant comme les limites 
d’autres intégrales qui n’ont pas cet inconvénient. C’est dans 
cette classe qu’on doit ranger celles dont M. Bidone a donné 
les valeurs dans les Mémoires de Turin pour l’année 1812(1). 

3 °. Lorsque fx renferme un radical, il devra conserver le 
même signe, ou plus généralement, avoir pour facteur la 
même racine de l’unité, dans toute l’étendue de l’intégration; 
et de cette manière , l’intégrale aura le même nombre de 
valeurs différentes, réelles ou imaginaires, dont le radical 
sera susceptible. Si le radical passe du réel a l’imaginaire, 
les signes dont il sera affecté dans ces deux périodes de va¬ 
leurs, n’auront pas de dépendance mutuelle, et la valeur de 
l’intégrale sera nécessairement ambiguë,c’est-à-dire, qu’après 
avoir donné un signe déterminé à la partie réelle, on pourra 
supposer indifféremment que la partie imaginaire soit multi¬ 
pliée par -h l/-1 ou par —: il est inutile d’insister sur 
cette circonstance qu’il suffit d’avoir indiquée. Lorsque le 
radical sera constamment réel entre les limites de l’intégra¬ 
tion, la nécessité d’un signe constant, sera une condition 
essentielle qui influera sur l’expression de l’intégrale définie. 
Pour en donner un exemple connu, considérons l’intégrale : 


W Voyez aussi sur ce point les Exercices de calcul intégral , tome II, 
page 125. 
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sin. xdx 

Yy i — 2 a cos. x 4- 



dans laquelle a est une constante positive ; et convenons de 
regarder le radical contenu sous le signe J, comme une quan¬ 
tité positive dans toute retendue de l’intégration. Il faudra 
prendre pour sa valeur, i 4- cl h. la limite et i —a 

ou a — i à la limite a?=o, selon qu’on aura a < i ou a > i. 
D’après cela , on trouve ces deux expressions différentes : 


/ 
J O 


sin .xdx 


1/7 


- 2 a cos. x 4- a * 


= 2, OU =- 


la première ayant lieu dans le cas de a < i, et la seconde 
dans le cas de a > i ; et si l’on différentie cette intégrale 
par rapport à a, on aura, dans les memes cas, cet autre 
exemple : 

/»* (a — coa.x)ùn-xdx _ ou = £ . 

I 1 a * 

J o (i— *acos.x-\-a 3 ) » 


Tl est à remarquer que si l’on fait a= i, les deux valeurs 
de la première intégrale sont égales, mais non pas celles 
de la seconde. Dans ce cas particulier, la dernière intégrale 
a pour valeurs zéro ou 2, selon qu’auparavant on regardait 
a comme plus petit, ou comme plus grand que 1 unité. Ce 
paradoxe tient à ce que les deux valeurs que 1 on détermine 
de cette manière, sont celles qui répondent à la différence 
1 —a infiniment petite, positive ou négative, et qu’à cette 
limite, un changement infiniment petit dans la valeur de 
1— a, suffît pour produire un changement brusque dans 
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celle de l’intégrale dont il est question. Dans le cas où l’ on 
aurait rigoureusement a = i,la valeur de cette intégrale serait 
la moyenne des deux valeurs précédentes, ou égale à l’unité; 
et, en effet, on a, dans cette hypothèse, 



(a — cos.ar )sin. xdx 

5 

i —a a cos. x -+- a*)~* 



Au reste, il existe beaucoup d’autres intégrales définies, 
renfermant, comme celle-ci, une constante sous le signe^ 

qui ont des valeurs différentes, selon que cette constante est 
positive ou négative, bien quelle puisse être infiniment pe¬ 
tite. C’est ainsi qu’on a, par exemple, 



selon que la constante a est > o, =o, ou < o. 


